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Abstract 

Almost nothing is known about the layer structure of solutions to strongly coupled 
systems of convection-diffusion equations in two dimensions. In some special cases 
we present first results. 
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1 Introduction 

In the survey paper [1] the authors present some results for weakly coupled systems in 
one and two space dimensions but state concerning strongly coupled convection-diffusion 
problems ,,we have only a limited grasp of the situation”. We aim to provide an inside 
into the nature of such problems at least in some special cases. 

A practical example of strongly coupled systems of convection-diffusion equations in 
2D (related to magnetohydrodynamic duct flow) is numerically studied in [2], namely 

^ -sAu + aVb = fi, 

—eAb -|- aS/u = f 2 

with some boundary conditions. 

Let us more generally consider the vector-valued function u = (mi,M 2 )^ solving the 
system 

c)lI ! 

(1.2a) -sAu + Ai- -hA 2 ^- \-pu = f in D, 

UXi C/X 2 

(1.2b) M = 0 on r = dVt, 

where e is a small positive parameter. We assume the matrices Ai,A 2 to be symmetric 
and and that the unit outer normal z/ = (z/i, 1 / 2 ) to D exists a.e. on dQ. 

Because 

E ^ ^ j{i^-Au, u)dr - f ((div A)u, u) 
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it is standard to assume 


^ dxi ’ 


(1.3) p >-sup II divAlloo- 

Then, problem fll.2p admits a unique weak solution. 

To describe the reduced problem we introduce the matrix 

(1.4) B := UiAi + z/ 2^2 

Suppose B to be nonsingular, i.e., dQ to be noncharacteristic. Then B allows the decom¬ 
position 

(1.5) B = B+ + B-, 

where B~^ is positive semidehnite, B~ negative semidehnite and the eigenvalues of B~^ are 
the positive eigenvalues of B and 0. The reduced problem to fll.21) is then given by 

(1.6) + puo = f with = 0 on T. 

OXi 0 X 2 

In [12] it was proved that u converges for e —)■ 0 to Mq- But concerning the convergence 
rate we only know the result of [5] for a problem with different boundary conditions: for 
/ G one has in the L 2 norm 

(1.7) ll«-Kollo <C£'/^II/II.. 

In the literature not much is known about the structure of layers. In the next Section we 
will discuss this question in, at least, some special cases and consider problems with two 
small parameters, as well. 

2 The reduced problem and location of layers 

Let us hrst revisit the situation in ID studied in [8]. We take the example 

(2.1) -3) “'=()) 

with m(0) = m( 1) = 0. In the one parameter case the eigenvalues and eigenvectors of A 
determine the asymptotic structure of the solution. Here, the eigenvalues are 5 and -5 with 
the corresponding eigenvectors (2,1)^ and (1, —2)^. With that knowledge and the general 
solution of the reduced equation we made in [S] the Ansatz 

(2.2) = wo + Wh + di exp(-5a;/e) -h ^2 (^^ 2 ) “ ^)/^) • 








The four constants ci,C 2 ,di,d 2 are computed from the boundary conditions. Eliminating 
all exponentially small terms one gets indirectly the solution of the reduced problem, for 
our example 


Al/25a;-8/25\ 
V-2/25 a; - 4/25yl ' 


Observe that Uq does not satisfy either of the given boundary conditions. 

There are two possibilities to derive the boundary conditions for the reduced equation 
directly. First as discussed in the Introduction, we can decompose the matrix B. 

In our example we have 


{Av, v) = ?,vl + %ViV 2 - 3^2 = (2ni + V 2 Y - (ui - 2^2)^- 


Therefore, the reduced problem is given by 




with the boundary conditions 

(2.3) (2Mj + n2)(i) = 0 and (uj - 2 m2)(0) = 0. 


An alternative approach consists in the diagonalization of the matrix A. Set 

u = Tv or V = T~^u. 


After transformation T ^AT is a diagonal matrix, we obtain a decoupled system in v. In 
our example, 

u = Tv with T = f ^ ^ ) 

\/5Vl -2; 

generates 

(2.4a) -ev'l - 5v[ = 

(2.4b) -ev2 + 5^2 = 

/^From the sign of the coefficients of the hrst order derivatives one can conclude where the 
layers of v are located. In our example vi has a layer at a; = 0, ^2 at x = 1. The back 
transformation yields: 

• the correct boundary conditions for the reduced problem are given by (12. 3 h 

• both components Ui and U 2 of the solution have layers at x = 0 and x = 1. 
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The transformation techniqne can also be applied for special systems in two dimen¬ 
sions. Let us assume that the constant matrices Ai,A 2 of the system fll.2p admit the 
representation 

(2.5) A^ = T diag(Al, A?)) T-\ A2 = T diag(A^, A^)) T"' 

with the orthogonal matrix 

sin (j) cos (j) 

— cos (j) sin 0 

for some (j) G [0, 27r]. Remark that this is the case if and only if ^1^42 = ^ 2 ^ 1 - 
Then, the transformation 



u = Tv or V = T 

yields also a decoupled system, namely (for p = 0, but this is not essential) 

(2.6a) + + = A 

(2.6b) + + 

That means: the sign of (A{, A^)^ • v determines the location of the layers of Ui, the sign 

of (A^, A^)^ ■ the location of the layers of V 2 . Let ns introduce 

r+ = {xer: (Aj,Ai) • j/> 0} 

and analogously Tj”, T)*, T^, T^, T^. Then, the conditions 

{T~i^ denotes the hrst row of T~^) are the correct boundary conditions for the reduced 
equation. 

Assume hi = (0,1)^. Then we have three typical cases (see Figure [I]). 

• (i) all A are positive 

Then ui and U 2 have overlapping layers at xi = 1 and at X 2 = 1. 

• (ii) only A^ is negative 

Then we have again overlapping layers at X 2 = 1, but ordinary layers of ui and U 2 at 
xi = 0 and at xi = 1. 

• (hi) both components of Ai and A 2 have a different sign 

Then ui and U2 have ordinary layers at every edge of the unit square. 
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Figure 1: Location of boundary layers 


Let us remark that the duct flow example mentioned in the introduction under the 
assumption a < 0 of [ 2 ] falls under case (hi) with 


T 


1 
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Let us now consider the two parameter case 

//?/ //?/ 

(2.7a) -SAu + Ai- -- \-pu = f 

uXi UX2 

(2.7b) u = 0 


in 

on r = dil 


with S = diag(ei,£ 2 ) and two small positive parameters. We are mostly interested in the 
case El /62 1 . 

Let us hrst collect some facts from [S] for the one-dimensional system 


(2.8) — diag(£i, e 2 )u" — Bu' + Au = /, m(0) = m(1) = 0. 


Assume 


B 



with bii > 0 and D = &11622 + ^12 > 0- 


Surprisingly, not the eigenvalues of B generate the structure of the boundary layers. The 
layer at a; = 0 is characterized by the exponential exp(—fen/^i), but the layer at x = 1 by 
the exponential exp(^A)_ 

The quantities —611 and are quotients of the leading principal minors of the matrix 
B. These values arise as well in the LDLA decomposition of the matrix B. In accordance 
with that observation the authors of [ 6 ] proposed to introduce new variables vhj v = L'^u 
resulting in 


— 8 *v" — diag((ii, d 2 )v' + A*v = /*, n( 0 ) = n(l) = 0 


* 


-L-^SiL-y. 


(2.9) 

with 


S 
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Unfortunately, the matrix 8 * is not diagonal. But its structure 


8* = 


El —El I 

-El I /^ + E2 


tells us (remember E 1 /E 2 1) the following: Assuming, for instance, di < 0 and ^2 > 0, 

both solution components have strong layers at x = 0 of the structure exp (((ii/£i)x). 
However, at x = 1 a weak layer of V 2 = M 2 of the structure 0 {ei/e2 ) exp ( — {d 2 /E 2 ){)- —x)) 
generates a strong layer of Vi and, consequently, Ui. In [9] this solution behavior was 
derived from some asymptotic approximation. 

The signs of di,d 2 also dehne the boundary conditions of the reduced problem, under 
our assumptions 


( 2 . 10 ) 

Consider the example 


■Uo(l) = 0, and Mo(0) +/mo(0) = 0- 


diag(£i,£2)M" 


u = 


0 


with m(0) = m( 1) = 0. Here 

^" (4/3 r 

and D = diag(3, —25/3). Consequently, in the case Ei E 2 we conclude 
• the boundary conditions for the reduced problem are 

Mg(l) = 0, and Mo(0) + 4/3 mq(0) = 0 


the solution u is characterized by the layer terms 

0(e 


^^exp(-^x) and exp ( - ^(1 - x)). 
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Remark 1 In the one parameter case the decomposition fll.5p allows to define the boundary 
conditions for the reduced problem. But with two parameters, for instance in ID, the matrix 
8 ~^B is not symmetric. Symmetrization of that matrix with a transformation of the type 


T-\£-'B)T with T = f) 


and adequately chosen fi is possible. But it is simpler to apply the LDL^ decomposition 
described above. □ 
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The results for the two parameter problem fl2.8l) can now be generalized to fl2.7p if we 
assume 


(2.11) Ai = L diag((i}, df) A 2 = L diag((i 2 , d^) ■ 

The transformation 


yields the system 


( 2 . 12 ) 


u = or V = u 

-£iAni - eiq/Av2 + d^ ■ Vvi + p{vi + av2) = /i, 
-Siq/Avi - [siq^ + £2)An2 + d^ ■ Vv2 + pv2 = h- 


Next we introduce as above rh,r^,r^,r 2 with respect to the vectors d^,d‘^. Both com¬ 
ponents Vi and V 2 as well as Ui and U 2 have strong exponential layers related to ei on T^. 
But on T^ the component V 2 = U 2 has only a weak layer. That weak layer generates a 
strong layer of ui related to the larger parameter 62 - 

The correct boundary conditions for the reduced equation are 

= 0 and n 2 |r- = 0- 

Summarizing: Under some restrictive conditions it is possible to characterize the possi¬ 
ble layers of strongly coupled convection-diffusion systems in 2D. Moreover, the adequate 
boundary conditions for the reduced problem can be derived (important for numerical 
methods which hrst solve the simpler reduced problem). 
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